Plant Molecular Biology

Plant regeneration is a major determinant of transformation and genome-editing efficiency, yet the endogenous regulatory networks controlling regenerative competence in horticultural crops remain incompletely understood. The miR319-TCP module regulates multiple developmental processes in plants, but its function in lettuce regeneration has not been defined. Here, we performed a genome-wide analysis of the TEOSINTE BRANCHED1/CYCLOIDEA/PROLIFERATING CELL FACTOR (TCP) gene family in lettuce (Lactuca sativa). Thirty-three LsTCP genes were identified and classified into Class I/PCF, Class II/CIN, and Class II/CYC/TB1 groups. Five CIN-class genes, LsTCP2, LsTCP3, LsTCP4, LsTCP10, and LsTCP24, were predicted as high-confidence miR319 targets and supported by degradome-based cleavage evidence. MIR319-overexpression (OX319) explants showed enhanced de novo shoot regeneration, with 94.5% regeneration efficiency and 1.92 shoots per explant, whereas STTM-miR319 suppression (S319) explants showed reduced regeneration, with 28.5% regeneration efficiency and 0.36 shoots per explant. These phenotypes were associated with altered expression of several miR319-targeted CIN-TCP genes, particularly LsTCP4, LsTCP10, and LsTCP24. Disruption of LsTCP4 increased regeneration efficiency to 91.4% and shoot production to 2.05 shoots per explant, resembling the regeneration-enhancing effect of miR319 overexpression. In contrast, disruption of the non-target CIN gene LsTCP17 did not significantly affect regeneration under the tested conditions. Together, these results identify LsTCP4 as a key miR319-responsive negative regulator of de novo shoot regeneration and highlight miR319-mediated repression of LsTCP4 as a potential endogenous strategy for improving lettuce regeneration.

The TETRATRICOPEPTIDE THIOREDOXIN-LIKE (TTL) proteins are a plant-specific family proposed to function as peripheral membrane proteins that contribute to abiotic stress tolerance in Arabidopsis, likely by maintaining cell wall integrity through brassinosteroid signaling. Previously, we identified a TTL gene that was differentially regulated at the translational level in nitrogen-fixing soybean plants under water deficit (WD) conditions. This finding prompted the characterization of the soybean TTL gene family. Using the Glycine max v4.0 proteome, we identified ten TTL homologs (GmTTL1-GmTTL10), which are unevenly distributed across five chromosomes. Phylogenetic and structural analyses grouped these genes into three clades and revealed a highly conserved exon-intron organization. Likewise, GmTTL proteins display a conserved number and arrangement of TPR and TRXL motifs. To gain insights into their potential biological functions, we integrated co-expression and differential expression analyses. This approach identified a co-expression module enriched for translationally downregulated genes related to the Gene Ontology terms "cellular anatomical entity", "membrane", "cell periphery", "cell wall modification", "nitrate assimilation", and "cell wall organization or biogenesis". Protein-protein interaction network analysis of this specific subset of genes uncovered a novel GmTTL connection with two nitrate reductase enzymes in nitrogen-fixing plants subjected to WD, potentially linking the TTL gene family to new functions or roles. This study provides a framework for future functional studies of GmTTL proteins and their contribution to abiotic stress adaptation in soybean. Key MessageThis work presents the first functional characterization of TTLs proteins in legume species and highlights key processes that may link the TTL gene family to new functions or roles.

Background and AimsPaspalum notatum reproduces through either sexuality or apomixis, two pathways that may coexist within the same individual and are regulated by interconnected molecular networks responsive to environmental cues. Here, we characterized the transcript structure and expression of BUD SITE SELECTION PROTEIN 13 (BUD13), a component of the RES spliceosomal complex previously reported as differentially expressed in florets of sexual and apomictic plants, as a first step toward testing its involvement in the molecular regulation of the apomixis-sexuality switch. MethodsPreviously generated floral and leaf transcriptomes from sexual and apomictic Paspalum notatum plants, including Oxford Nanopore long-read data, were mined to characterize BUD13 transcript structure and expression. Phylogenetic analyses and in silico mapping were conducted to infer evolutionary relationships and determine the origin of the transcripts. Differential expression was validated by RT-qPCR, while in situ hybridization was used to reveal cell-specific ovule expression patterns. Key resultsBUD13 is expressed in Paspalum notatum florets as a truncated isoform (SHORT) encoding a small protein lacking part of the herpes simplex virus regulatory protein (ICP4) domain. Two SHORT transcripts, SHORT1 and SHORT2, with different 5' untranslated region (UTR) regions, were identified in flowers. SHORT1 was consistently upregulated in apomictic ovules from premeiosis to anthesis. Both transcripts originated from a single genomic locus located in the subtelomeric region of the short arm of chromosome 6. SHORT isoforms with variable structures were detected in other monocots. In situ hybridization showed that, whereas BUD13 was expressed throughout sexual ovules, expression was absent from the female germline of apomictic ovules. A consistent expression was observed in somatic proembryos of aposporous embryo sacs. ConclusionsOur findings reveal structural, spatial and temporal divergence in BUD13 expression between sexual and apomictic reproductive programs, providing new insights into the molecular regulation of asexual seed formation.

Ovate Family Proteins (OFPs) are a class of plant-specific, negative nuclear transcriptional regulators characterized by conserved C-terminal OVATE domain. This study on comparative functional characterization of two head-to-head arranged OFPs - AtOFP2 (Ovate-OFP with full ovate domain) and AtOFP17 (Ovate-Like OFP with partial ovate domain) provides critical insight into how structural variations in ovate domain leads to functional divergence. Detailed phenotypic analysis of 28 physical and physiological traits of loss- and gain-of-function mutants revealed that both genes act as broad, pleotropic repressors of plant growth and development. Removal of repression in knock-down mutants of both genes exhibited reduced duration of seed dormancy, faster rate of germination and growth, bigger plants and significantly higher seed yield. In contrast, constitutive over-expression showed a generalized repressive nature of both genes, with nuanced differences for fine tuning of specific traits. For example, both genes showed antagonistic behaviours on root hair architecture. AtOFP2 act as a strong repressor of root hair development whereas AtOFP17 is a stronger repressor of hypocotyl and root cell architecture. AtOFP17 owing to partial ovate domain exerts a mild level of repression throughout life span as indicated by smaller plants and lesser yield in knock-down AtOFP17 mutants. On the contrary, AtOFP2 exerted a much stronger repressor effect in which > 90% over-expression mutants died at the juvenile stage ; the survival of remaining 10% is probably owing to activation of dosage-dependent feedback loop mechanism as indicated by normal growth of mature plants, and is also evident by transcriptome data. Transcriptome analysis of roots of 7-day old seedling of knock-down and over-expression mutants of AtOFP2 showed downregulation of OFP2 in over-expressed mutants. However, severely stunted phenotype indicated presence of stable OFP2 protein to exert effects. Analysis of DEGs in OFP2 mutants revealed that it acts as an important regulator working at intersection of hormonal signalling affecting critical genes required for auxin, cytokinin, GA, BR and ABA functioning. Perturbations across hormonal signalling pathways affects cell wall remodelling factors such as EXPANSINS, Xyloglucan hydrolases (XTHs) and cellulose synthases (CSLs) causing overall stunted growth; and epidermal patterning genes such as WER, GL1, EGL3, TTG1 leading to severely reduced root length and root hairs. Significantly, functional analysis of this master regulator highlighted a significant economic potential. Knockdown of both these genes relieves their natural repression on reproductive traits, leading to longer siliques, bigger and heavier seeds, and substantially increased overall seed yield, positioning AtOFP2 and AtOFP17 as highly valuable targets for agricultural crop improvement.

O_LIEngineering nitrogen fixing crops requires not only transferring the nitrogenase structural genes, but also the accessory genes to synthesize its iron-sulphur cofactors. Scaffold protein NifU is a critical element in this system as the starting point of nitrogenase cofactor assembly. NifU has been successfully produced in plants, however, its optimal production required high levels of iron in the medium. This is likely due to a faulty connection with the endogenous iron trafficking network C_LIO_LITo identify specific elements targeting iron to NifU, pull-down assays were performed to identify showing bacterioferritin A (BfrA) as a likely candidate. Co-immunopurification, mutant characterization, iron transfer assays, and co-expression in Nicotiana benthamiana assays were carried out. C_LIO_LIBfrA transfers iron to NifU through protein-protein interactions. When these two proteins were co-expressed in N. benthamiana leaves, there was an increase in NifU production. In turn, it led to doubling NifH synthesis, a nitrogenase structural protein that is also required for the synthesis of the more complex nitrogenase cofactors. C_LIO_LIOur results provide a new element towards engineering nitrogen-fixing crops. They also underscore the importance of transferring the metal delivery systems when expressing metalloproteins in heterologous systems. C_LI

Potato leafroll virus (PLRV) is an economically important viral disease of potato (S. tuberosum). Genetic resistance to this phloem-limited virus is rare, and no cloned resistance (R) genes have been reported. Rladg confers resistance to PLRV in an Andean potato landrace, LOP-868 (Velasquez et al. 2007). We identified the functional Rladg gene as a homolog of the tomato TIR-NLR-encoding Bs4. Rladg interacts with the serine protease domain of the PLRV protein P1, which is essential for virus replication. This recognition is independent of the proteases enzymatic activity, and the Rladg immune receptor oligomerizes upon direct association with the protease. Like PLRV, many poleroviruses contain a serine protease. Despite their diverse amino acid sequences, these proteases are predicted to share similar structures. Rladg recognizes all ten tested polerovirus proteases, suggesting a conserved structural recognition mechanism. We propose that Rladgs broad recognition capacity could enable resistance to poleroviruses in many crop species. Rladg is the first R-gene reported to confer resistance to a phloem-limited pathogen and could provide enhanced resistance to many economically important poleroviruses.

Elongating rice internodes present a developmental gradient from dividing meristem to mature cells, providing an elegant pseudo-time course for study of plant vegetative development. We tested the hypothesis that DNA methylation regulates gene expression during rice internode development by integrating RNA-seq and bisulfite DNA sequencing across eight internode segments. Previously described topologically associated chromatin domain borders aligned with transcription start sites of constitutive expressed genes. CpG and CHG differential methylation was enriched in young segments, consistent with maintenance methylation; whereas CHH methylation showed similar differential abundance in young and old segments. CHH and CHG methylation in upstream regions, CpG methylation within gene bodies, and any methylation in 5' and 3' untranslated regions were permissive of moderate to high gene expression. Very low expression was associated with CpG methylation upstream, CHG and CHH methylation within gene bodies, and CpG and CHG methylation downstream. A nonrandom subset of genes, including cell wall-related glycoside hydrolases, lignin and tricin biosynthesis enzymes, and WD40 proteins, showed methylation-expression correlations, with expression changes enriched in triple-marked elements. These results suggest that internode phenotypes of DNA methylation machinery mutants relate to alteration of specific target genes, opening approaches for grass culm improvement for lodging resistance and biomass production.

In Asteraceae species, the capitulum is a compact inflorescence, featuring a characteristic reproductive structure. Despite the identification of a few key regulatory factors, the transcriptome-level information on the developing capitulum remains limited. Here, we applied single-cell and spatial transcriptome sequencing to investigate the developing Gerbera hybridas capitulum during floret differentiation. We obtained a transcriptomics atlas encompassing different stages of the Gerbera capitulum and analyzed the cellular and spatial dynamics of gene expression. Using marker gene expression and GO enrichment of cluster-specific DEGs, we annotated putative cell types and described changes in gene expression across sampled stages, potentially associated with ongoing developmental processes. We detected activity of previously undescribed MADS-box genes and defined their spatial expression patterns. Notably, the MADS-box gene GAGL12 was found to be enriched in the putative capitulum phloem cells. The GAGL12 protein was shown in yeast two-hybrid assays to interact with several other MADS-domain proteins with hypothesized functions in vasculature development, and further detailed in silico analyses supported a candidate role in the development of capitulum vasculature. Altogether, we provide integrative and dynamic transcriptomic insight into capitulum and floret development and lay a basis for future functional studies of the control and development of this intriguing reproductive structure.

Transgenerational plasticity can shape plant phenotype and influence plant response to environmental changes in interaction with the current conditions. While how past stress interact with either current optimal or stress conditions is increasingly documented within a single plant, transgenerational plasticity remains particularly poorly understood especially at the metabolome level. In our study, we investigated whether heat stress induces transgenerational metabolic and phenotypic modifications along two successive clonal ramet generations of the sub-Antarctic aquatic plant Limosella australis. We performed untargeted metabolomic approaches and measured morphologic and performance traits, to assess both transgenerational plasticity of the metabolome and the phenotype. We found that heat stress remodelled the metabolic profile and influenced the foraging strategy of our clonal plant, and that some of these metabolic changes persisted into the first clonal generation. This one therefore adopted an intermediate growth strategy, even though culture conditions were optimal. By comparing differentially accumulated features between daughter ramets from heat stressed mother ramets and from unstressed mother ramets, we identified common and specific metabolites accumulation to heat stress response, belonging to diverse compound families. However, we did not observe any adaptative advantage and any metabolic imprint during another heat stress applied on the second clonal generation. This work provides especially new clues into how plant metabolome integrates and transfers previous stressed clonal generation's information.

The receptor kinase FERONIA (FER) perceives endogenous RAPID ALKALINIZATION FACTOR (RALF) peptides and regulates a plethora of plant physiological processes, including immunity. RALF peptides also bind to LEUCINE-RICH REPEAT EXTENSIN (LRX) proteins as structural components of the cell wall. We recently showed that the FER-RALF-LRX pathway supports colonization by the obligate biotrophic fungal pathogen Erysiphe cruciferarum (Ecr), a member of the powdery mildew species complex that infects Arabidopsis. Genetic disruption of the pathway primarily affects conidiation of the fungus, raising the question of effects on fungal nutrition. To get further insight into the underlying mechanisms, we performed RNA sequencing (RNAseq) to identify differential transcriptional responses of FER-RALF-LRX pathway mutants upon Ecr infection. Surprisingly, our results revealed that pathway disruption has a limited impact on the overall transcriptional changes upon fungal infection. However, consistent with previous reports, FER-RALF-LRX pathway mutants show changes in basal expression of a plethora of genes, mainly associated with cell wall metabolism, jasmonic acid signalling, amino acid biosynthesis and secondary metabolism. Many of these genes are regulated by Ecr infection across genotypes, too. This raises the question whether these are relevant pathway components for powdery mildew host establishment downstream of the FER-RALF-LRX module. In summary, our data reveals new insights into FER-RALF-LRX-dependent responses that may support host susceptibility to biotrophic plant pathogens.

Rice (Oryza sativa) is highly sensitive to salinity, yet the metabolic mechanisms underlying salt tolerance remains incompletely understood. In this study, we performed leaf tissue-specific untargeted metabolomic profiling of the salt-tolerant introgression line JN100 (JN), its donor parent Nona Bokra (NB), and its recurrent parent Jupiter (JU) to characterize metabolic responses to salt stress. Comparative analysis identified differentially accumulated metabolites (DAMs) spanning diverse chemical classes, including amino acids, sugars and carbohydrates, lipids, organic acids, cofactors, electron carriers, and nucleotides. Under salt stress (SS), 201 DAMs (89 upregulated and 112 downregulated) were detected in JN relative to JU. Notably, metabolites such as allantoin, glycitin, nicotinamide ribotide, D-arabinono-1,4-lactone, violanthin, L-methionine S-oxide, ribitol, lysine, rutin, glutamine, pantothenic acid, and quinic acid, showed significant differential accumulation. Pathway enrichment analysis revealed significant enrichment of arginine biosynthesis, purine metabolism, and alanine, aspartate, and glutamate metabolism, indicating extensive reprogramming of nitrogen and energy-associated metabolic pathways under salinity stress. Integration of transcriptomic and metabolomic datasets from the SS experiments further identified ten differentially expressed genes (DEGs) associated with the metabolite network in the JN vs. JU comparison. Among these, OsDHQDT/SDH, OsFd-GOGAT, phenylalanyl-tRNA synthetase, OsP5CS1, OsP5CS2, and a pyridoxal phosphate-dependent transferase were linked to metabolites involved in shikimate, amino acid, and proline metabolism. Collectively, these results demonstrate that salinity tolerance in rice is associated with coordinated transcriptional and metabolic reprogramming that supports oxidative stress mitigation and adaptive stress responses.

Crop photoperiodism models assume that flowering time is primarily controlled by daylength, yet many field observations contradict this view. We previously proposed an alternative framework integrating daily changes in sunrise and sunset times (dSR and dSS). Variety trials in Madagascar and in Argentina supported this concept: mid-late sorghum varieties from the northern hemisphere flowered late or very late when sown in November and December, consistent with the higher dSR/dSS values of the southern hemisphere summer. One Malian variety, sown monthly over six years in West Africa, exhibited high interannual variability in flowering time when sown between November and February. This revealed that up to four photoperiodic responses -- two quantitative and two qualitative, occurring at different times of the year -- may coexist within a single late photoperiod sensitive variety. All responses use only dSR and dSS cues. The qualitative responses are triggered by an internal phasic coincidence, which is set by a linear relationship between dSR and dSS at the onset of plant photoperiod sensitivity, and between dSR+dSS at panicle initiation. The research model fitted data from 28 varieties grown in Mali well. It also accurately fitted the duration to PI observed in three varieties sown at tropical and temperate latitudes. HighlightThe seasonal photoperiodic adaptation of flowering time in sorghum plants may rely on several signal transduction pathways regulated by sunrise and sunset times rather than day length.

Developmental processes depend on carefully coordinated gene expression. Expression is modulated by the binding of transcription factors (TFs) to cis-regulatory elements (CREs), like enhancers and promoters. Many computational and experimental approaches have been developed to find CREs, particularly enhancers, in the genome, each with strengths and caveats. Given the increasing availability of ATAC-seq data and methods to find TF binding therein, we hypothesized that we could use TF footprinting tools to find clusters of TF binding events within accessible chromatin that may act as CREs. Using Drosophila anterior-posterior patterning network as a test bed, we used a digital genomic footprinting tool (DGT), TOBIAS, on previously published early embryo ATAC-seq data to characterize the TF footprint landscape of 16 TFs essential for embryonic patterning. Even in this system, with its extensive enhancer annotation, most footprinted TF binding sites lie outside of known enhancers, with intergenic and intronic regions hosting the highest TF footprint count, albeit at low density. To find potential novel enhancers, we identified high-density TF footprint clusters that are highly conserved and overlap with active enhancer histone mark signals. Five high confidence candidates were selected for reporter assay validation and all five were found to drive spatially patterned expression in the embryo. This study shows that even in a highly characterized system, the analysis of footprinted TF binding sites in ATAC-seq data can uncover new regulatory regions and suggests this approach may be helpful in using existing ATAC-seq data to find novel CREs. ARTICLE SUMMARYGiven the increasing availability of ATAC-seq datasets, workflows to exploit the data to uncover new cis-regulatory elements (CREs), including enhancers, are valuable. Using early anterior-posterior patterning in the Drosophila embryo as a test case, we find that previously published transcription factor footprinting tools and ATAC-seq data can be analyzed to yield new candidate CREs. Experimental validation confirms the activity of selected candidate CREs, suggesting that existing data can be analyzed to find novel regulatory elements.

Codon usage bias is a fundamental genomic characteristic that prefers non-random preferential use of synonymous codons. It is a major determinant of translational efficiency, gene regulation, and molecular evolution. However, the evolutionary bias and functional relevance of codon usage bias across the plant lineage is poorly defined and yet to understand what are the major factors responsible for relative synonymous codon usage (RSCU) in genomes and how codon usage bias influences the gene regulation, molecular evolution genomes. A genome-wide codon usage bias study of coding DNA sequences of 262 plant genome was conducted. It encompassed more than 4.6 billion codons from > 11 million coding sequences. Relative synonymous codon usage, codon adaptation index, codon-anticodon mapping, effective number of codon (ENC)-GC3, GC1,2-GC3, parity rule 2 (PR2-bias), molecular economy, and machine learning approaches were used for the study. It was found that codon usage bias was strongly non-random and exhibited a clear phylogenetic structuring. The higher plants favoured A/T-ending, whereas early-diverging lineages were enriched in G/C-ending codons. Analysis of RSCU, codon adaptation index, and codon-anticodon pairing indicated that translational selection is mediated by tRNA availability, contributing sustainability to these molecular patterns. Machine-learning approaches identified a small subset of codons having outsized influence on genome-wide codon usage landscapes. Further studies revealed the presence of robust inverse relationships between the effective number of codons and GC content at synonymous third positions. Neutrality analysis revealed approximately 61% of variation was driven by mutational pressure, tempered by selective constraints. Phylogenetic reconstruction showed a progressive relaxation of codon bias from algae to angiosperms while maintaining a conserved molecular economy cost of ~ 30 ATP per codon across the lineages. The study revealed codon usage bias is lineage-specific evolutionary conserved trait governed by mutation, selection, and translational optimization.

Translation initiation signals shape gene expression across all domains of life. In eukaryotes, nucleotide constraints surrounding the start codon are commonly described by the Kozak Consensus Sequence (KCS), whereas in bacteria and archaea, initiation frequently involves Shine--Dalgarno ribosome-binding motifs. Although these signals have been extensively characterized in model organisms, their large-scale diversity and evolutionary distribution remain incompletely explored. We present KozakExplorer, a reproducible framework for quantitative and comparative analysis of translation initiation contexts from genome assemblies and annotations. The software performs strand-aware extraction of start codon environments from FASTA and GFF3 files and applies information-theoretic metrics---including Kullback--Leibler (KL) divergence and information content (IC)---to measure positional nucleotide constraints relative to a background model. Derived summary statistics (Kozak Strength Index [KSI], maximum information content, peak position) convert motif patterns into interpretable per-genome signatures suitable for cross-species comparison. Our primary analysis covers 2,282 eukaryotic reference genomes, producing a standardized dataset of translation initiation metrics. Dimensionality reduction via t-SNE on per-position KL divergence, information content, and motif nucleotide frequencies reveals a structured eukaryotic KCS landscape with kingdom-level clustering and continuous variation in signal strength. A dedicated case study of 216 Apicomplexa genomes shows genus-level structure consistent with host range and phylogeny. An extended analysis across 25,344 reference genomes (22,253 bacteria, 809 archaea) places eukaryotic patterns in a global comparative framework, revealing transitions between sharply localized Kozak motifs and distributed Shine--Dalgarno-type signatures. Implemented within the open-source Sequana ecosystem, KozakExplorer is distributed as a Python module and an interactive web application that accepts local annotated assemblies, GenBank records, or NCBI RefSeq accessions, and exports all computed metrics, embeddings, and coordinates for downstream comparative and evolutionary genomics.

Auxin is a key phytohormone that regulates all aspects of plant growth, development, and environmental responses, making the precise analysis of its distribution and signaling essential for understanding plant adaptation and physiological processes. However, despite the agricultural importance of oilseed rape (Brassica napus), the lack of robust, species-specific molecular tools limits detailed studies of hormone signaling in this crop. Here, we developed and characterized reporter systems for the sensitive visualization and quantification of auxin distribution and signaling in B. napus. The DR5cc auxin signaling reporter and a novel synthetic auxin-responsive reporter, BIP3, assembled from promoter fragments of three oilseed rape IAA genes, were generated to drive GUS expression. In hairy roots, both reporters showed auxin-responsive expression in the root apical meristem that became broader after auxin treatment. In transgenic seedlings, flowers at anthesis, and 12-day-old embryos, DR5cc exhibited a more defined expression pattern than BIP3. To monitor real-time auxin dynamics under abiotic stress, DR5cc fluorescent reporters were employed in hairy roots. Mannitol and NaCl treatments induced a time-dependent increase in fluorescence, peaking at 6-12 h before returning to basal levels after 24 h. Furthermore, dual-reporter assays enabled simultaneous monitoring of auxin and cytokinin signaling, revealing distinct hormone-specific spatial responses in hairy roots. Finally, we established a quantitative DII (qDII) reporter system using degron domains from B. napus Aux/IAA proteins, providing a high-resolution quantitative readout of auxin depletion. Together, these reporter systems enable spatial, temporal, and quantitative analyses of auxin dynamics during development and stress adaptation in oilseed rape.

Lima bean (Phaseolus lunatus) is a broadly adapted, economically important leguminous crop and a susceptible host of root-knot nematodes (Meloidogyne spp.; RKN), which are a devastating plant pathogen in agricultural systems worldwide. To date, there have been few studies to elucidate the genetic determinants of RKN resistance in lima beans. Understanding the genetic mechanisms underlying resistance is essential for improving resistance traits and incorporating them into lima bean breeding programs. To assist in marker-assisted selection, we aimed to identify and map quantitative trait loci (QTLs) conferring RKN resistance-related traits. Three recombinant inbred line (RIL) populations were used in this study. Three populations were derived by crossing two RKN-resistant parents with the same RKN-susceptible parent and with each other. All populations were genotyped using genome-wide single-nucleotide polymorphism (SNP) markers. Each population was screened for root galling (RG) and RKN egg reproduction (ER) in response to M. incognita and M. javanica in greenhouse experiments. Three major QTLs were detected and mapped on chromosome Pl04 (QRk-pl04.1), Pl05 (QRk-pl05.1) and Pl10 (QRk-pl10.1) across populations. Among them, QRk-pl05.1 and QRk-pl10.1 affected levels of RG and ER of both RKN species, while QRk-pl04.1 suppressed root galling and reproduction responses of M. incognita but not of M. javanica. These chromosomal regions defined by flanking markers will help guide marker-assisted breeding and gene discovery for broad-based RKN resistance in lima beans.

Plant regeneration is a dynamic developmental process that spans from cell dedifferentiation to organ reconstruction in response to inductive cues, such as wounding stress and hormonal signals. Although this capacity varies widely both between and within species, a comprehensive understanding of the genetic and epigenetic bases of this variation remains incomplete. To address this issue, we revisited published datasets on natural variation in in vitro regeneration capacity in Arabidopsis thaliana. Using quantitative genetic approaches, including meta-analyses of genome-wide association studies (GWAS) and multi-locus models, we dissected the genetic architecture underlying regeneration traits. Our results showed that shoot regeneration capacity is primarily explained by allelic variation in the cis-regulatory region of WUSCHEL (WUS), a key regulator of shoot meristem formation. Notably, these polymorphisms are also associated with epigenetic variants of the DNA transposon ATDNA2T9C, which is located within the regulatory region. Furthermore, allelic variation in ARABIDOPSIS RESPONSE REGULATOR 2 (ARR2), a positive regulator of cytokinin signaling, is associated with callus formation and greening traits and may promote shoot formation through genetic interactions with WUS alleles. Although in vitro regeneration is controlled by complex, multilayered gene regulatory networks, our results suggest that, in A. thaliana, natural variation in regeneration capacity is largely shaped by a small number of major-effect modifiers together with epigenetic variation and genetic interactions, despite the substantial heterogeneity observed among natural populations.

Heavy metals act as essential metalloprotein cofactors in numerous physiological processes but can become toxic when non-essential metals accumulate or when essential metals are in excess. As plants continuously encounter heavy metals through their roots, they have evolved complex homeostatic mechanisms to regulate metal uptake and distribution. The Heavy Metal ATPase (HMA) gene family encodes a group of heavy metal transporting P-type ATPases that have been linked to stress resistance and nutrient supply. Here, we used a bioinformatics approach to identify and characterise 13 HMA genes containing characteristic P1B-type ATPase domains and motifs in the barley Morex V3 reference genome. The genes are located on five of the seven barley chromosomes. Phylogenetic analysis revealed that they cluster into five sub-clades, including one clade unique to barley. Expression profiling across multiple datasets showed distinct temporal and tissue-specific expression patterns among HvHMAs, with several members exhibiting significant transcriptional responses to specific biotic and abiotic stresses. By utilising recently available pan-transcriptomic and pan-genomic resources, we have identified substantial allelic diversity and inter-accession variation in HvHMAs. Our findings suggest that HvHMAs have functions extending beyond canonical heavy metal homeostasis and warrant further investigation for their potential roles in broader physiological and stress-related processes.

Plant PHAS locus discovery remains challenging because phasiRNA-producing loci must be distinguished from other sRNA-producing regions with high abundance or apparent periodicity. This problem is especially acute for reproductive 24-PHAS loci, which occur within genomes that also produce abundant 24-nt siRNAs from non-PHAS regions. We present Phasis, an open-source Python software tool for plant PHAS-locus discovery from small RNA sequencing data. Phasis combines statistical evidence for phased accumulation with locus-level features and a Register-Resolved Locus Interpretation Layer that evaluates whether candidate loci show coherent phased architecture. Across diverse plant datasets, Phasis recovered validated or annotated 21- and 24-PHAS loci with a strong balance between call-level precision and reference-locus recall, and generally outperformed PhaseTank and ShortStack in matched benchmark analyses. The register-resolved interpretation layer reduced unsupported calls by separating coherent phased loci from ambiguous sRNA-producing regions. In maize dcl5 mutant libraries, Phasis showed strong depletion of 24-PHAS recovery, supporting DCL5-dependent recovery of reproductive 24-PHAS signal. Together, these results support Phasis as a biologically interpretable tool for large-scale discovery of plant DCL-dependent phasiRNA loci.